New aspects of microwave properties of Nb in the mixed state 
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Abstract 

We present a study of the frequency dependence of the vortex dynamics in a conventional superconductor. We have 
employed a swept- frequency, Corbino-disk technique to investigate the temperature (S.GK-Tc) and high-field (from 
Hcijl to HcTz) microwave complex resistivity in Nb thin (20-40 nm) films as a function of the frequency (1-20 GHz). 
We have found several previously unnoticed features: (i) a field-dependent depinning frequency in the GHz range; (ii) 
deviations from the accepted frequency dependence [1], that can be ascribed to some kind of vortex creep; (iii) the 
presence of switching phenomena, reminiscent of vortex instabilities. We discuss the possible origin of the features 
here reported. 
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1. Introduction 

The microwave response of superconductors in the 
mixed state is a particularly suitable probe to in- 
vestigate the short-range vortex dynamics. Histor- 
ically, vortices in conventional superconductors are 
thought to follow the well-established Gittleman- 
Rosenblum predictions [T], based on single- vortex 
response to an alternating current without thermal 
effects. Surprisingly, there are only a few experi- 
mental reports on this behaviour, in particular as 
a function of frequency, temperature and magnetic 
field. The advancement on knowledge about high 
frequency vortex dynamics in the recent years [2], 
triggered by the discovery of high- Tc superconduc- 
tors, and the revamped interest in conventional su- 
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perconductors as potential components of devices 
based on superconductor/ferromagnet multilayers 
pi, stimulates a direct, careful experimental deter- 
mination of the vortex dynamics in conventional su- 
perconductors. 

When vortices are set in motion by time- varying 
currents, they experience heterogeneous forces 
which include a viscous drag force, phenomenologi- 
cal representation of vortex motion-induced power 
dissipation, pinning forces, arising from interaction 
with material defects and hindering vortex move- 
ment, and stochastic forces of thermal origin which 
promote vortex detachment from pinning sites. 
The electrodynamic high frequency response aris- 
ing from the interaction of the fluxon system with 
microwave currents has been modeled by many au- 
thors |l|2|4|5j . On very general grounds, all models 
can be represented through a universal expression 
for the complex vortex resistivity [6] : 
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Pvm — Pvm,l H~ ^Pvm,2 



where the resistivity pff represents the free flux 
flow value, reached in the high frequency limit in 
which vortices experiences only the dissipative vis- 
cous drag, vq is a characteristic frequency and the 
dimensionless parameter < e < 1 is a measure of 
the weight of creep (depinning events) phenomena. 
We adopt the Coffey Clem (CC) model [2], in which 
e = [/o(w)]"^ and i^o = i^pil - e)~^ (/i(u)//o(w)). 
Here, u = Ua{T, B)/{2KbT) is the normalized en- 
ergy barrier height within the assumption of a uni- 
form periodic pinning potential of height C/q, Vp is 
the so-called pinning frequency, which marks the 
crossover between the vortex pinned motion and 
the purely dissipative motion, arising at ^ fp, 
/„ is the n-order modified Bessel functions of the 
first kind. The CC model reverts to the simpler 
Gittleman and Roscmblum (GR) expression = 

Pff (1 + i (vp/v)) / (l + i^p/'^)^^ M for vanishing 
creep (e —> 0), so that vq — > Vp. In the latter, the 
pinning frequency can be immediately derived as 

Hvm,2/ yvm,l • 

The exact determination of the three vortex pa- 
rameters comes directly from comparison to the 
measured frequency dependence, while only esti- 
mates can be obtained from single-frequency mea- 
surements resorting to an accurate analysis [6 . 
Neglecting the weight of some of the parameters 
in order to simplify the analysis may significantly 
affect the estimates of the parameters. 

In the following, we illustrate microwave measure- 
ments performed on Nb thin films by means of the 
Corbino disk technique, which allows wide band, 
frequency dependent determinations of the complex 
resistivity. 



2. Experimental results and discussion 

We measure the swept frequency {^=1-20 GHz) 
microwave response of Nb samples, at low and con- 
stant input power (<0.1 mW). Temperature is var- 
ied in the range 3.6 K<T <Tc (stability < 0.001 K) 
through a Helium flow cryostat. A static magnetic 
fleld Hc2/2 < H < Hc2 is applied perpendicularly 
to the sample surface. 

The studied Nb thin films have thicknesses d—2Q- 
40 nm and were grown on 0.5 mm thick, 5x5 mm^ 
square AI2O3 substrates in a ultra high vacuum 
dc diode magnetron sputtering system (10^^ nibar 



base pressure; lO^'^ mbar sputtering Argon pres- 
sure) . The fabrication was realized at room tempera- 
ture. The deposition rate of 0.3 nm/s was controlled 
by a quartz crystal monitor calibrated by low-angle 
X-Ray reflectivity measurements performed using 
a Philips X-Pert MRD high resolution diffractome- 
ter. The sample thickness was determined by fitting 
the refiectivity profile of the sample with a simula- 
tion curve obtained following the Parrat and Nevot- 
Croce formalism [718] . The fit revealed the forma- 
tion of a Nb205 oxide layer of the order of 2-3 nm 
on the top of film which causes a reduction of the 
effective Nb layer thickness. This in turn results in a 
depression of Tc, which is strongly thickness depen- 
dent for low thickness values such as ours 9j. Four 
contact measurements of the dc resistivity on a typi- 
cal sample 20 nm thick, similar to the one here stud- 
ied, yielded a constant normal state p„=22 /xficm 
for temperatures Tc < T < 30 K. The critical tem- 
perature Tco = 6.37 K, evaluated as the midpoint of 
the resistive transition. 

The microwave response is measured through 
a Corbino disk [iD] setup: a swept frequency mi- 
crowave radiation, generated by a Vector Network 
Analyser (VNA) system, is fed in a coaxial cable 
short-circuited on the sample. We use a launcher 
with a custom-made, spring loaded central pin in 
order to ensure good contact with the Nb film. 

We measured the (complex) reflection coefficient 
r(i') of the electromagnetic wave impinging upon 
the sample surface. The modulus |r|, bounded 
within [0,1], is a good indication of power dissipa- 
tion of the (super)conducting sample. In particular, 
a vanishingly low dissipation (as in the Meissner 
state) would yield |r| = 1, whereas by increasing 
dissipation |r| would decrease below 1. A quantita- 
tive analysis requires the extraction of the complex 
effective surface impedance Z^ffiv) of the sample 
according to the well-known expression [llj : 



et)\ ) "1 -r(z^) 



(2) 



where Zq is the characteristic impedance of the coax- 
ial line. The small thickness of our samples allows the 
use of the thin film approximation in which Zg/ / = 
p/d, where p is the complex resistivity of the su- 
perconductor. We did not observe, in the [T, H , v] 
region here explored, resonances such as the well 
known substrate resonances occurring in thin films 

The actually measured reflection coefficient 
T„i^{v)^ determined by the VNA at its input, include 



2 



the sample response T(i>) as well as the interposed 
coaxial line response. The calibration of the line 
contribution is a critical issue, as described in Ref. 
[lOj . Here, we used as reference the normal state so 
that all results are presented as normalized resistiv- 
ities p/pn- 

In the following wc consider a Nb sample 20 nm 
thick (Nb20). Raw data are presented in Fig. [T] to 
demonstrate the effect of a magnetic field and of the 
measuring frequency on the superconducting tran- 
sition. The black circles correspond to single fre- 
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Fig. 1. |r| vs T at selected fields and frequencies, sample 
Nb20. Circles: = 5 GHz; squares: U2 = 17 GHz. 

quency data taken at vi = b GHz and poH = T 
at various T. A finite width of the transition is due 
to finite quasiparticle contribution that can be ob- 
served at high frequency. 

The application of a magnetic field (white cir- 
cles in Fig. [1]) determines an evident and progres- 
sive broadening of the transition, signature of the 
increased dissipation caused by the motion of the 
vortices penetrated in the sample. 

In a magnetic field, dissipation increases at larger 
frequency (y2 = 17 GHz, squares in Fig.[T]), directly 
showing that the vortex characteristic frequency lies 
in our measuring frequency range. 

If, as customarily performed, we rely on single fre- 
quency measurements and derive the characteristic 
frequency neglecting creep (GR model, Ref. 
we obtain different values by using data at different 
frequencies: 2.1 GHz using data at and 3.8 GHz 



^ Vp can be computed directly from the complex F, 
given the straightforward derivation Pvrn,2{'^)) / Pvm.,1 = 
2|F| sin((/))/{l - |F|2) where </> = arg(r). 



using data at 1^2 ■ Clearly, such a procedure is unsuit- 
able, and wide band measurements are required. 

A sample measurement in the full 1-20 GHz range 
is reported in Fig. [51 taken at T = 3.66 K and poH = 
1.3 T. One first, striking, feature is the presence 
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Fig. 2. Normalized p\ and Ap2 (black and grey symbols, 
respectively) vs v, sample Nb20. Thin continuous line: fit 
according to CC model (see text). 

of switching-like phenomena appearing simultane- 
ously on both the real, pi, and imaginary, Ap2 = 
P2(H) — p2{Q), parts of the complex resistivity. Most 
likely the phenomenon originates from vortex insta- 
bilities, as we will discuss in the following. Here we 
only stress that it can be thought as superimposing 
to an underlying, unperturbed response, represented 
by the lower/upper envelopes of the real and imag- 
inary parts, respectively, of the measured j5. The 
measured (unperturbed) pi + iAp2 shows no signs 
of quasi-particle/superfluid pairbreaking contribu- 
tions (as can be expected not too close to Bc2{T)): 
pi increases and then saturates at high v whereas 
Ap2 presents a maximum. These features are clear 
signatures of single vortex dynamics, so that it can 
be safely taken pi -\- iAp2 = Pvm- 

Therefore, a fit of pi +iAp2 can be done by follow- 
ing the CC model. Fit results are reported as thin 
continuous lines in Fig. [2| their quality is remark- 
able. The obtained fitting parameters are Pff/pn = 
0.7, Vp = 7.2 GHz, e = 0.41. 

It can be seen that the pinning frequency Vp dif- 
fers significantly from the simple GR estimate previ- 
ously calculated. This fact is connected to the pres- 
ence of significant creep processes {e = 0.41 upon 
a theoretical maximum value of 1), which evidently 
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excludes the use of simple, GR-like models. This is 
an important result of this work: we found that sig- 
nificant creep is present even at low (3.66 K) temper- 
atures. It is worth stressing again that these results 
are made possible by the use of wide band measure- 
ments, which demonstrates to be essential for the 
comprehensive study of vortex dynamics. 

Moreover, a reliable determination of vortex pa- 
rameters enables also any successive comment on 
their eventual field dependence. By increasing the 
field, we observe an increase of creep and a decrease 
of Up^ which can be explained by a weakening of vor- 
tex pinning by approaching the Hc2{T) line. 

Finally, the obtained pff yields, within the 
Bardeen Stephen model (p// = PnB/Bc2) [E], a 
value of Bc2 1.9 T compatible with the rough 
estimate of i?c2 ~ 2.2 T (taking B k, iiqH), deter- 
mined as the field at which the superconducting 
signal vanishes within the sensitivity of our system. 

We now come back to the switching phenomenon. 
This is the second main result of this work. In- 
deed, to the best of our knowledge, this is the first 
time that similar phenomena are observed in the 
microwave range. Our whole set of measurements 
shows that this phenomenon appears for H > Hc2 / 2 
and increases in amplitude as H Hc2- Moreover, 
switches occur upward on pi and downward for p2'. 
this feature points to commutations of the vortex 
system forth to and back from a higher dissipation 
state. In order to get additional insight, we per- 
formed continuous wave, fixed frequencies measure- 
ments. A sample of them is reported in Fig. [3l in 
terms of (uncalibrated) \ measured as a function 
of time in a sample = 30 nm thick (Tc ~ 7.5 K). The 




0.376 




t(s)10 11 12 14 16 18 20 22 24 



Fig. 3. ITml vs time in residual field iiqH ~ 30 mT, Nb 
sample 30 nm thick. 

time-dependent nature and the pseudo-periodicity 
of the switches is clearly apparent, as well as the 
"large" time scales involved (pseudo-period=1.2-1.5 



s) . Continuous wave data also show no apparent de- 
pendence on both the frequency and the power of 
the microwave stimulus. Finally, the phenomenon 
has been observed in all the samples and resulted 
independent from the instrumentation used. 

Although the origin of this switching is at the mo- 
ment unclear, two main factors are relevant. First, 
the noted correlation to the vicinity of the Hc2{T) 
line suggests some form of thermally-driven, fluc- 
tuation enhanced metastability involving the vor- 
tex system [15]. Second, the Corbino disk geometry 
likely plays a role, since the strong velocity gradi- 
ents it induces on the fluxon lattice are considered 
the source of peculiar and unconventional vortex dy- 
namics regimes |16I17| . 

3. Summary 

We have studied the frequency dependence of the 
vortex dynamics in Nb thin fllms at different tem- 
peratures and magnetic fields in the range 1-20 GHz 
through the Corbino disk technique. The resulting 
single vortex dynamics was well modelized within 
frameworks such the Coffey Clem one. In particular, 
all the three main vortex parameters were clearly 
determined, showing the presence of significant 
creep contribution even at our lowest (3.66 K) tem- 
peratures. In addition we observed time-dependent 
switching phenomena likely related to metastabili- 
ties involving the fluxon system. This phenomenon 
is probably linked to thermal and fluctuation effects 
involving the vortex system, as well as to the pecu- 
liar fluxon velocity pattern imposed by the Corbino 
disk geometry. Further studies will be necessary to 
ascertain the exact nature of the phenomenon. 

This work has been partially supported by an Ital- 
ian MIUR-PRIN 2007 project. 
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